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Abstract 

The neutralino-nucleon cross section in the context of the MSSM with universal 
soft supersymmetry-breaking terms is compared with the hmits from dark 
matter detectors. Our analysis is focussed on the stability of the corresponding 
cross sections with respect to variations of the initial scale for the running of 
the soft terms, finding that the smaller the scale is, the larger the cross sections 
become. For example, by taking 10^°~^^ GeV rather than Mqut, which is a 
more sensible choice, in particular in the context of some superstring models, 
we find extensive regions in the parameter space with cross sections in the range 
of 10^^-10^^ pb, i.e. where current dark matter experiments are sensitive. For 
instance, this can be obtained for tan/3 > 3. 
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1 Introduction 



Recently there has been some theoretical activity [|T], 0, ^, ||, |], |^ [] analyzing the com- 
patibility of regions in the parameter space of the minimal supersymmetric standard 
model (MSSM) with the sensitivity of current dark matter detectors, DAMA |[ll|] and 
CDMS [0]. These detectors are sensitive to a neutralino-nucleon cross section^ (T-j^o_p 
in the range of 10~^-10~^ pb. Working in the supergravity framework for the MSSM 
with universal soft terms, it was pointed out in |^, Q that the large tan/5 regime 
allows regions where the above mentioned range of o'^o_p is reached. Besides, working 
with non-universal soft scalar masses, they also found <y^o^_p 10^^ pb for small val- 
ues of tan/5. In particular, this was obtained for tan/3 > 25 (tan/3 > 4) working with 
universal (non- universal) soft terms in The case of non- universal gaugino masses 
was also analyzed in with interesting results. 

The above analyses were performed assuming universality (and non-universality) 
of the soft breaking terms at the unification scale, Mqut ~ 10^^ GeV, as it is usually 
done in the MSSM literature. Such a scale can be obtained in a natural manner within 
superstring theories. This is e.g. the case of type I superstring theory [|^, [l^ and 
heterotic M-theory [0, |T^. However, recently, going away from perturbative vacua, 
it was realized that the string scale may be anywhere between the weak scale and the 
Planck scale. For instance D-brane configurations where the standard model lives, 
allow these possibilities in type I strings [|1^, [1^, |I9|, Similar results can also 



be obtained in type II strings pTI and weakly and strongly coupled heterotic strings 



22| , |23| . Hence, it is natural to wonder how much the standard neutralino-nucleon 
cross section analysis will get modified by taking a scale Mj smaller than Mqut for 
the initial scale of the soft terms 0. 

The content of the article is as follows. In Section 2 we will briefly review several 
scenarios suggested by superstring theory where the initial scale for the running of the 
soft terms is Mj instead of Mqut- In particular, we will see that Mi ^ \{]^'^-^'^ GeV is 
an attractive possibility. The issue of gauge coupling unification, which is important for 
our calculation, will also be discussed. Then, in Section 3, we will study in detail the 



iSec also @, |, |, |To) for other works. 

^Let us recall that the lightest neutralino x\ is the natural candidate for dark matter in supersym- 
metric theories with conserved R parity, since it is usually the lightest supersymmetric particle (LSP) 
and therefore stable 

•^This question was recently pointed out in for a different type of analysis. In particular 
the authors studied low-energy implications, like sparticle spectra and charge and colour breaking 
constraints, of a string theory with a scale of order 10^^ GeV. Similar phcnomenological analyses were 
carried out in the past pq| for Mpianck rather than AIqut- 
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stability of the neutralino- nucleoli cross section with respect to variations of Mj. For 
the sake of generality, we will allow Mj to vary between 10^^ GeV, which corresponds 
to Mguti and 10^" GeV. Of course the results will be valid not only for low-scale string 
scenarios but also for any scenario with an unification scale smaller than Mqut- Let us 
finally remark that the analysis will be carried out for the case of universal soft terms. 
This is the most simple situation in the framework of the MSSM and can be obtained 
e.g. in superstring models with dilaton-dominated supersymmetry breaking or in 
weakly and strongly coupled heterotic models with one Kahler modulus Finally, 
the conclusions are left for Section 4. 



2 Initial scales 

As mentioned in the Introduction, it was recently realized that the string scale is not 
necessarily close to the Planck scale but can be as low as the electroweak scale. In this 
context, two scenarios are specially attractive in order to attack the hierarchy problem 
of unified theories: a non-supersymmetric scenario with the string scale of order a few 
TeV [|18|, [19|, and a supersymmetric scenario with the string scale of order 10^'^^^^ GeV 
H. Since we are interested in the analysis of supersymmetric dark matter, we will 
concentrate on the latter. 

In supergravity models supersymmetry can be spontaneously broken in a hidden 
sector of the theory and the gravitino mass, which sets the overall scale of the soft 
terms, is given by: 

F 

~ — , (1) 

J-VJ- Planck 

where F is the auxiliary field whose vacuum expectation value breaks supersymmetry. 
Since in supergravity one would expect F ^ Mpi^^^i^, one obtains m3/2 ~ Mpianck 
and therefore the hierarchy problem solved in principle by supersymmetry would be 
re-introduced, unless non-perturbative effects such as gaugino condensation produce 
F ^ My/Mpianck- Howcver, if the scale of the fundamental theory is M/ 10^''"^^ 
GeV instead of Mpianck, then F ^ Mj and one gets m3/2 ~ M^y in a natural way, 
without invoking any hierarchically suppressed non-perturbative effect . 



For example, embedding the standard model inside D3-branes in type I strings, the 
string scale is given by: 

Mf = '^Mfl^Ml , (2) 
V2 

where a is the gauge coupling and is the compactification scale. Thus one gets 
Mi ^ 10i°-i2 GeV with M^ ^ 10^-^° GeV. 
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There are other arguments in favour of scenarios with initial scales M/ smaller than 
Mgut- For example in |22[ scales Mi ^ lO^o-i^ GeV were suggested to explam many 



experimental observations as neutrino masses or the scale for axion physics. These 
scales might also explain the observed ultra-high energy 10^° eV) cosmic rays as 
products of long-lived massive string mode decays. Besides, several models of chaotic 
inflation favour also these initial scales p8| . 



Inspired by these scenarios we will allow the initial scale M/ for the running of the 
soft terms to vary between 10^^ GeV and 10^'^ GeV, when computing the neutralino- 
nucleon cross section below. As we will see, the values of the gauge coupling constants 
at those scales will be crucial in the computation. This head us for a brief discussion 
of gauge coupling unification in models with low initial scale: 

(a) Non universality of gauge couplings 
An interesting proposal in the context of type I string models was studied in [15, 29|: 



if the standard model comes from the same collection of D-branes, stringy corrections 
might change the boundary conditions at the string scale M/ to mimic the effect of 
field theoretical logarithmic running. Thus the gauge couplings will be non universal 
and their values will depend on the initial scale M/ chosen. This is schematically 
shown in Fig. la for the scale M/ = 10^^ GeV, where g^ ^ 0.8, g2 ~ 0.6 and gi ~ 0.5. 
Clearly, another possibility giving rise to a similar result might arise when the gauge 
groups came from different types of D-branes. Since different D-branes have associated 
different couphngs, this would imply the non universality of the gauge couplings. 

(b) Universality of gauge couplings 
On the other hand, if gauge coupling unification at Mj, = a, is what we want 
to obtain, then the addition of extra fields in the massless spectrum can achieve this 
task ^0[. An example of additional particles which can produce the beta functions. 



63 = —3, 62 = 3, 61 = 19, yielding unification at around M/ = 10^^ GeV was given in 



2 X [(1, 2, 1/2) + (1,2, -1/2)] + 3 X [(1, 1,1) + (1, 1, -1)] , (3) 

where the fields transform under S'[/(3)c x SU {2)l x U{1)y- In this example one obtains 
g{Mi) ^ 0.8. This is schematically shown in Fig. lb. 

As a matter of fact, once the Pandora's box containing extra matter fields is open, 
other possibilities arise. Note that the example above does not contain extra triplets, 
and therefore runs as in the MSSM 

' ' i^iln^. (4) 



a3(Q) a:,{M,„,„) 2ir M 



susy 
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where 63 = —3 and Msusy indicates the supersymmetric threshold. However, introduc- 
ing extra triplets, 63 = —3 + ^^3 will increase and therefore 0:3 will also increase. 
Likewise, extra doublets and/or singlets will allow to increase the value of 02 and ai 
as in the example above and therefore we will be able to obtain unification at Mj, but 
for bigger values of = a. For example, the additional particles 

3 X [(3, 1, 2/3) + (3, 1, -2/3)] + 6 x [(1, 2, 1/2) + (1, 2, -1/2)] , (5) 

produce 63 = 0, 62 = 7, bi = 27, yielding unification again at around M/ = 10^^ GeV 
but for g{Mj) ^ 1.3. 

We will see in the next section that due to the different values of the gauge couplings 
at Mj, scenarios (a) and (b) give rise to qualitatively different results for cross sections. 



3 Neutralino-nucleon cross sections versus initial 
scales 

In this section we will consider the whole parameter space of the MSSM with the 
only assumption of universality. In particular, the requirement of correct electroweak 
breaking leave us with four independent parameters (modulo the sign of the Higgs 
mixing parameter /i which appears in the superpotential W = fiHiH2). These may 
be chosen as follows: m, M, A and tan (3, i.e. the scalar and gaugino masses, the 
coefficient of the trilinear terms, and the ratio of Higgs vacuum expectation values 

<H2> 

<Hi>- 

On the other hand, we will work with the usual formulas for the elastic scattering 
of relic LSPs on protons and neutrons that can be found in the literature 0, ||, ||. 



In particular, we will follow the re-evaluation of the rates carried out in 0, using their 
central values for the hadronic matrix elements. 

As mentioned in the introduction, the initial boundary conditions for the running 
MSSM soft terms are usually understood at a scale Mqut- Smaller initial scales, 
as for example Mj ^ 10^^ GeV, will imply larger neutralino-nucleon cross sections. 
Although we will enter in more details later on, basically this can be understood from 
the variation in the value of /^^ with Mi, since the cross sections are very sensitive to 
this value. 

Let us discuss then first the variation of /i^ with M/. Recalling that this value is 
determined by the electroweak breaking conditions as 

,2 _ mj,^-ml^tan^(3 1 . 
^ " tan2/3-l 2 ^' 
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we observe that, for tan/5 fixed, the smaller the initial scale for the running is, the 
smaller the numerator in the first piece of becomes. To understand this qualitatively, 
let us consider e.g. the evolution of m|^^ (neglecting for simplicity the bottom and tau 
Yukawa couplings) 



where t = \n{Mf/Q^ 



mjj^it) = + G{t) , (7) 



and the functions fi{t) are given by: 

/ 



1 72 I • (9) 



V 



An bit 



For Mj = Mgut "we recover the usual case of the MSSM: t 61 for Q ~ 1 TeV. 
However, for Mj smaller than Mgut-, the value of t will decrease and therefore fi{t) 
will also decrease, producing a smaller value of G{t). As a consequence m^^ at Q ~ 1 
TeV will also be smaller. For rn'\j^{t) the argument is similar: the smaller the initial 
scale for the running is, the less important the negative contribution m?^^ to /x^ in (H) 
becomes. 

In fact, when the initial scale is decreased, it is worth noticing that the values of 
the gauge couplings are modified. Therefore, this effect will also contribute to modify 
the soft Higgs mass-squared (see e.g. (|^)). Let us consider first the ) with non- 

universal gauge couplings at M/ discussed in the previous section (see also Fig. la), 
where Q!2(M/) and ai{Mi) are smaller than a{MGUT)- For instance, for this has 
obvious implications: fi(t) decrease for M/ smaller than Mgut, not only because t 
decreases but also because ai{0) are smaller. 

The above conclusions have important consequences for cross sections. As is well 
known, when \fi\ decreases the Higgsino components, Ni^ and of the lightest 
neutralino 

X? = NuB"^ + iVi2iy° + NnHi + Ni^H^ (10) 

increase and therefore the spin independent cross section also increases. We show both 
facts in Figs. 2 and 3. 

These figures correspond to /i < 0. Opposite values of /i imply smaller cross sections 
and, moreover, well known experimental constraints as those coming from the 6 — 57 
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process highly reduce the /i > parameter space. The other parameters are chosen 
as follows. For a given tan/3 and neutralino mass, the common gaugino mass M is 
essentially fixed. For the common scalar mass m we have taken m = 150 GeV. Finally, 
for the common coefficient of the trilinear terms A we have taken A = —M. This 
relation is particularly interesting since it arises naturally in several string models 



26| , p!5| . In any case we have checked that the cross sections and our main conclusions 
are not very sensitive to the specific values of A and m in a wide range. In particular 
this is so for | A/M \ < 5 andQ 50 GeV < m < 250 GeV. 

We have checked that our results are consistent with present bounds coming from 
accelerators and astrophysics. The former are LEP and Tevatron bounds on supersym- 
metric masses and CLEO b s'y branching ratio measurements. The latter are relic 
neutralino density bounds and will be discussed in some detail later on. 

In Fig. 2, for tan/3 = 10, we exhibit the gaugino-Higgsino components-squared 
N?, of the LSP as a function of its mass m^o for two different values of the initial 
scale, Mj = 10^^ GeV ^ Mqut and Mi = 10^^ GeV. Clearly, the smaller the scale 
is, the larger the Higgsino components become. In particular, for Mj = 10^^ GeV 
the LSP is mainly Bino since A'^n is extremely large. The scattering channels through 
Higgs exchange are suppressed (recall that the Higgs-neutralino-neutralino couplings 
are proportional to N13 and A^i4) and therefore the cross sections are small as we will see 
explicitly below. As a matter of fact, the scattering channels through squark exchange 
are also suppressed by the mass of the first-family squarks. Indeed in this limit the 
cross section can be approximated as 




where rrir is the reduced mass and rriq, 6 are the mass and the mixing angle of the first- 
family squarks respectively. However for Mj = 10^^ GeV the Higgsino contributions 
A^i3 and A^i4 become important and even dominant for m^o < 130 GeV (e.g. with 
tan/5 = 3 this is obtained for m^o < 65 GeV). Following the above arguments this will 
imply larger cross sections. Indeed scattering channels through Higgs exchange are now 
important and their contributions to cross sections can be schematically approximated 
as 

^x?-P « T^^I^H (^7'iVii - 92^2) P , (12) 
^1 An m- 



h 



^Let us remark however that for m in the range 50 — 100 GeV the neutraUno is not the LSP in the 
whole parameter space. In some regions the stau is the LSP. 
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where i = 3,4, Xq are the quark Yukawa couphngs and mh represent the Higgs masses. 
It is also worth noticing that, for any fixed value of Mj, the larger tan/3 is, the larger 
the Higgsino contributions become. The reason being that the top(bottom) Yukawa 
coupling decreases (increases) since it is proportional to This implies that 

the negative (positive) contribution m'jj^^m'jjj to /x^ is less important. The discussion 
of the cases with tan /3 > 10 is more subtle and will be carried out below. 

The consequence of these results on the cross section is shown in Fig. 3, where the 
cross section as a function of the LSP mass m^o is plotted for five different values of 
the initial scale Mr. For instance, when m^o = 100 GeV, o"^o„„ for Mr = 10^^ GeV is 

' Ai ' A.1 y 

two orders of magnitude larger than for Mqut- In particular, for tan /? = 3, one finds 
cr^o_p < 10"'^ pb if the initial scale is M/ = 10^*^ GeV. However cr^o_p > 10"*^ GeV is 
possible if M/ decreases. For M/ < 10^^ GeV, taking into account the experimental 
lower bound on the lightest chargino mass m^± = 90, the range 70 GeV < m^o < 100 
GeV is even consistent with the DAMA limits. As discussed above, the larger tan (3 is, 
the larger the Higgsino contributions become, and therefore the cross section increases. 
For tan /5 = 10 we see in Fig. 3 that the range 60 GeV < m^o < 130 GeV is now 
consistent with DAMA limits. This corresponds to M/ < 10^^ GeV. 

Finally, we show in Fig. 3 the case tan/5 = 20. Then the above range increases 50 
GeV < m^o < 170 GeV, corresponding now to M/ < 10^^ GeV. It is worth noticing 
here that the value of /i^ is very stable with respect to variations of tan/5 when this 
is large (tan/? > 10). This is due to the fact that /i^ ~ ~™H2 ~ (^^e (§)). Since 
sin/5 ^ 1, the top Yukawa coupling is stable and therefore the same conclusion is 
obtained for Tn\j^ and /i^. Thus, for a given Mi, the reason for the cross section to 
increase when tan (3 increases cannot be now the increment of the Higgsino components 
of the LSP. Nevertheless there is a second effect in the cross section which is now the 
dominant one: the contribution of the down-type quark Yukawa couplings (see (p!^)) 
which are proportional to 

In Fig. 4 we plot the neutralino-proton cross section as a function of tan (3. Whereas 
large values of tan/? are needed in the case M/ = Mqut to obtain cross sections in 
the relevant region of DAMA experiment, the opposite situation occurs in the case 
Mi = 10^^ GeV since smaller values are favoured. 

Let us consider now the case (b) with gauge coupling unification at Mj discussed 
in Section 2 (see also Fig. lb). The result for the cross section as a function of m^o is 
plotted in Fig. 5 for tan (3 = 20. Clearly, the cross section increases when M/ decreases. 
However, this increment is less important than in the previous case. The reason being 
that now 0:2 (M/) and ai{Mi) are bigger than a{MGUT) instead of smaller. For example 
this counteracts the increment of fi{t) in (^ due to the smaller value of t when Mj 
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is smaller. Due to this effect, only with tan (3 > 20 we obtain regions consistent with 
DAMA limits. 

The results for the case with extra triplets (see e.g. (H)) are worst since the gauge 
coupling at the unification point M/ is bigger than above. The value of CTj^o-p at 



may be even smaller than its usual value at Mqut- 

Before concluding let us discuss briefly the effect of relic neutralino density bounds 
on cross sections. The most robust evidence for the existence of dark matter comes 
from relatively small scales. Lower limits inferred from the fiat rotation curves of spiral 
galaxies |3l|] are Vthaio ^ 10 Vt^s or Vthaio h"^ ^ 0.01 — 0.05, where h is the reduced 
Hubble constant. On the opposite side, observations at large scales, (6 — 20) Mpc, 
have provided estimates of QcDAih"^ ~ 0.1 — 0.6 but values as low as VLcDuh'^ ~ 
0.02 have also been quoted ||33|. Taking up-to-date limits on the baryon density from 
nucleosynthesis and overall matter-balance analysis one is able to obtain a favoured 
range, 0.01 < VLcDuh'^ ^ 0.3 (at ~ 2a CL) |3^, |3^. Note that conservative lower limits 
in the small and large scales are of the same order of magnitude. 

In this work, the expected neutralino cosmological relic density has been computed 
according to well known techniques (see ||13[)- In principle, from its general behaviour 
fi^/i^ oc l/{crann), whcre aann is the cross section for annihilation of neutralinos, it is 
expected that such high neutralino-proton cross sections as those presented above will 
then correspond to relatively low relic neutralino densities. However, our results show 
that for some of the largest cross-sections, with e.g. Mj ^ 10^^~^^ GeV, the value 
of the relic density is still inside the conservative ranges we considered above. As a 
function of the intermediate scale our results show a steady increase in the value of 
the relic density when we move from Mj ^ 10^^ GeV down to Mj ^ 10^^ GeV. In this 
respect, the main conclusion to be drawn from our results is that it should be always 
feasible, for large areas of supersymmetric parameters, to find an "intermediate" initial 
scale Ml which represents a compromise between a high neutralino cross section and 
an adequate relic neutralino density. 

We expect that neutralino coannhilitationsQ do not play an important role here 
since the mass differences among the LSP and the NLSP are not too small in general. 
For example, for the relative mass difference among the two lightest neutralinos, we 
have found Am^o^o/m^o > 0.2 in most of the interesting regions, being typically much 
higher. 



^The computation of the relic density carried out here, following p^ , contains only a partial 
treatment of neutralino coannhilation channels. 
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4 Conclusions 



In this paper we have analyzed the relevant implications for dark matter analyses of the 
possible existence of an initial scale Mj smaller than Mqut to implement the boundary 
conditions. 

We have noted that the neutralino-nucleon cross section in the MSSM is quite 
sensitive to the value of the initial scale for the running of the soft breaking terms. 
The smaller the scale is, the larger the cross section becomes. In particular, by taking 
Ml 10^°~^^ GeV rather than Mqut ~ 10^^ GeV for the initial scale, which is a 
more sensible choice e.g. in the context of some superstring models, we find that the 
cross section increases substantially being compatible with the sensitivity of current 
dark matter experiments cr^o_p ^ 10~^-10~^ pb, for tan/5 > 3. For larger values of 
the initial scale, as e.g. Mj = 10^^ GeV, the compatibility is obtained for tan/? > 10. 
Let us remark that these results have been obtained assuming non-universal gauge 
couplings at M/, as discussed in Section 2. They should be compared with those of the 
MSSM with initial scale Mqut, where tan /? > 20 is needed. 

We have also discussed the corresponding relic neutralino densities and checked 
that they are of the right order of magnitude in large areas of the parameter space for 
the neutralino being a CDM candidate. 

The above computations have been carried out for the case of universal soft terms. 
This is not only the most simple possibility in the framework of the MSSM, but also 
is allowed in the context of superstring models. This is e.g. the case of the dilaton- 
dominated supersymmetry breaking scenarios or weakly and strongly coupled heterotic 
models with one Kahler modulus. In this sense the analysis of neutralino-nucleon cross 
sections of those models is included in our computation. 

Obviously, non universality of the soft terms will introduce more flexibility in the 
computation, in particular in the value of /x^, in order to obtain regions in the parameter 
space giving rise to cross sections compatible with the sensitivity of current detectors. 
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Figure 1: Running of the gauge couplings with energy assuming non universahty (a) 
and universahty (b) of couphngs at the initial scale Mi — 10^^ GeV. Dashed hues 
indicate the usual running of the MSSM couplings. 
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Figure 2: Gaugino-Higgsino components-squared Nf.^ of the lightest neutrahno as a 
function of its mass, for the case tan j3 = 10, and for two different values of the initial 
scale. Mi = 10^^ GeV and Mj = 10^^ GeV. 



14 




10 
10 
10 
10 
10 



-10 



10 
10 
10 
10 
10 
10 
10 



-5 
-6 
-7 
-8 
-9 
-10 



10 
10 
10 
10 
10 
10 
10 



-10 



tanfS=5 



'immmimmMmimMimiiimiiii 




200 



m^^iGeV) 





tan/S=10 




DAMA 








~-=-:^r^7^ CDMS 




m,.=90GeV^\^ 





50 



100 



150 



200 



m^Q(GeV) 



tan/? = 20 
MMMMMMMMM. 



CDMS 



m„=90 GeV 



50 



100 



150 



200 



Figure 3: Neutralino-proton cross section as a function of the neutralino mass for 
several possible values of the initial scale Mj, and for different values of tan (3. Current 
experimental limits, DAMA and CDMS, are shown. The region on the left of the line 



denoted by m- 



90 GeV is excluded because of the experimental lower bound on 



the lightest chargino mass. 
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Figure 4: Neutralino-proton cross section as a function of tan (3 for two possible values 
of the initial scale, M/ = 10^^ GeV, M/ = 10^^ GeV, and for different values of the 
neutralino mass, namely 50, 75 and 100 GeV. DAMA limits are also shown. 
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Figure 5: The same as in Fig. 3 but only for tan/3 = 20 and with gauge couphng 
unification at M/. 
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